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We report on the mean-field study of the Chiral Magnetic Effect (CME) in static magnetic fields 
within a simple model of a parity-breaking Weyl semimetal given by the lattice Wilson-Dirac Hamil¬ 
tonian with constant chiral chemical potential. We consider both the mean-field renormalization of 
the model parameters and nontrivial corrections to the CME originating from re-summed ladder 
diagrams with arbitrary number of loops. We find that on-site repulsive interactions affect the chiral 
magnetic conductivity almost exclusively through the enhancement of the renormalized chiral chem¬ 
ical potential. Our results suggest that nontrivial corrections to the chiral magnetic conductivity 
due to inter-fermion interactions are not relevant in practice, since they only become important 
when the CME response is strongly suppressed by the large gap in the energy spectrum. 
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1. INTRODUCTION 

Anomaly-driven quantum transport phenomena in 
Dirac and Weyl semimetals, such as the Chiral Magnetic 
Effect (CME) [1] and the anomalous Hall effect [2], are 
nowadays in the focus of active theoretical and experi¬ 
mental studies. On the one hand, these phenomena are 
very interesting from a theorist’s point of view since they 
are nontrivial macroscopic manifestations of the universal 
laws of quantum physics at the microscopic level, much 
like the phenomena of superconductivity and superfluid¬ 
ity. On the other hand, implementation of anomalous 
transport in real materials would offer an intriguing pos¬ 
sibility of a dissipation-free transport of spin-polarized 
current. An experimental signature of the Chiral Mag¬ 
netic effect is the quadratic enhancement of the electric 
conductivity in the direction of the applied constant mag¬ 
netic field, which was originally predicted in [3] and has 
recently been observed in the Dirac semimetals Bii_a;Sba, 
[4], ZrTes [5] and NaaBi [6]. Very similar behavior of the 
electric conductivity in the direction of the magnetic field 
has also been found in holographic models [7, 8] and nu¬ 
merically in lattice QCD [9]. 

The Chiral Magnetic Effect is characterized by the chi¬ 
ral magnetic conductivity acME, which relates the in¬ 
duced electric current j and the external magnetic field 
B: 


j = CFCME B, 


( 1 ) 


where we assume the natural units for all the dimension¬ 
ful quantities. For free Dirac fermions, the chiral mag- 
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netic conductivity is [1] 


crcME = 


NfRA 

27r2 ’ 


( 2 ) 


where Nf is the number of Dirac cones in the energy 
spectrum and fj-A is the chiral chemical potential which 
parameterizes the imbalance of densities of left- and 
right-handed fermions. In experimental studies of the 
magneto-resistivity of Dirac/Weyl semimetals, this chi¬ 
rality imbalance is also generated due to the quantum 
anomaly and within a simple relaxation time approxima¬ 
tion can be estimated as [5] 


3avpT E ■ B 


( 3 ) 


where E is the external electric field which probes the DC 
conductivity of the material, is the Fermi velocity, 
T is the relaxation time of axial charge, and T is the 
temperature. Combining this formula with (1) and (2), 
it is straightforward to see that the electric conductivity 
should grow as a square of the magnetic field. 

The relations (1) and (2) can be also reproduced in 
various models of parity-breaking Weyl semimetals [10- 
15]. Strictly speaking, for lattice models this formula is 
only valid in the limit of a spatially uniform magnetic 
field which slowly changes in time [16] - in the case of 
a static spatially uniform magnetic field <jcme vanishes 
[17], in accordance with the general requirements of the 
gauge invariance of the effective action [18, 19]. How¬ 
ever, if the static external magnetic field is allowed to be 
spatially non-uniform, (Jcme becomes some non-trivial 
function of the wave vector k of the external field. More¬ 
over, if the chiral symmetry is not spontaneously broken, 
the asymptotic value of acME at large k is related 
to the axial anomaly coefficient in a nontrivial way [20] 
and approaches the universal value (2), but with a mi¬ 
nus sign. This change of sign as compared to the most 
commonly quoted result (2) is related to the vector cur¬ 
rent conserving regularization of (Jcme [20, 21] and does 





2 


not affect the relation between acME and the anomaly 
coefficient. 

So far, most calculations of the chiral magnetic con¬ 
ductivity were performed without taking into account the 
inter-electron interactions. This approximation is not un¬ 
reasonable, since in the hydrodynamics, kinetic theory or 
Fermi liquid approximations anomalous transport coeffi¬ 
cients take universal values even in the presence of inter¬ 
fermion interactions [22-26]. However, in recent years it 
has been realized that the chiral magnetic conductivity 
(as well as the other anomalous transport coefficients) 
can still get nontrivial corrections due to inter-fermion 
interactions in two important situations: first, when the 
chiral symmetry is spontaneously broken [20, 27] and sec¬ 
ond, when the electric current is coupled to a dynami¬ 
cal gauge field [27-30]. In the first case, the appearance 
of massless Goldstone modes violates the applicability 
of the hydrodynamical approximation, and the effective 
mass term which mixes chiralities makes the notion of a 
Fermi surface at finite ill-defined [20]. In the second 
case, the correlators of electric current j receive perturba¬ 
tive QED corrections, which result in a renormalization 
of acME- 

Clearly, both of these cases are relevant for real Weyl 
semimetals. Electrostatic inter-electron interactions in 
Dirac and Weyl semimetals are effectively enhanced by 
a factor of inverse Eermi velocity vp ~ 10“^ ... 10“^ 
[4, 5], thus the effective QED coupling constant Oe// 
oiqedIvf ~ 137 VF order of one. Therefore in 

Weyl semimetals one typically deals with a strongly cou¬ 
pled fermionic gas which might become unstable towards 
spontaneous chiral symmetry breaking. On the other 
hand, inter-electron interactions in the context of con¬ 
densed matter physics are always mediated by dynamical 
photon fields, therefore even interactions which are too 
weak to trigger spontaneous symmetry breaking can in 
general result in perturbative corrections to ucme- 

One should also mention that the classical dynamics of 
an electromagnetic field coupled to the plasma of chiral 
fermions at finite chiral chemical potential is likely 
to be unstable with respect to the decay of the chiral 
chemical potential at the expense of increasing magnetic 
helicity, as the studies within the frameworks of the chi¬ 
ral kinetic theory [31, 32], magnetostatics [33, 34] and 
the AdS/CET correspondence [35-37] suggest. There¬ 
fore modelling the parity-breaking Weyl semimetals with 
a time-independent Dirac-like Hamiltonian with finite ^a 
is a certain idealization valid at timescales much shorter 
than the typical decay time of chirality imbalance. The 
origin of this chiral plasma instability is the coupling to 
the magnetic field, which is suppressed as l/vp in Dirac 
and Weyl semimetals. Hence one can expect that the 
typical timescale of the decay of fiA is rather long, much 
longer than the timescales at which e.g. the formation of 
the chiral condensate due to strong electrostatic interac¬ 
tions takes place. Recent real-time simulations of chiral 
plasma instability [38] support this expectation. Another 
situation in which constant fj,A might be a reasonable ap¬ 


proximation is when the decay of chirality imbalance due 
to chiral plasma instabilities or other chirality-changing 
processes is compensated by a constant inflow of chirality 
in parallel electric and magnetic fields [5, 39]. Since in 
experiments one typically deals with a strong constant 
magnetic field [4-6] and an electric field of rather low 
frequency, a reasonably good approximation is to assume 
that the magnetic field and the electric current induced 
by it are also static. In this case, the calculation of the 
current becomes considerably simpler and can be carried 
out entirely in the equilibrium-state formalism. Another 
motivation for considering the static limit is that at suf¬ 
ficiently small frequencies and momenta the difference 
between the static and time-dependent chiral magnetic 
conductivities vanishes once one takes into account re¬ 
laxation processes [40]. 

In this paper we study the effect of strong inter¬ 
electron interactions on the static chiral magnetic con¬ 
ductivity in a parity-breaking Weyl semimetal. In par¬ 
ticular, we are interested in the question whether strong 
inter-electron interactions can lead to significant devia¬ 
tions of the chiral magnetic conductivity from the uni¬ 
versal value (2). A particularly interesting possibility is 
that the interactions can strongly enhance the CME re¬ 
sponse, as suggested by the mean-field analysis for the 
continuum Dirac fermions [27]. Understanding the influ¬ 
ence of inter-electron interactions on the chiral magnetic 
effect is thus important for the correct interpretation of 
the recent experiments [4-6]. 

However, the continuum Dirac Hamiltonian is only the 
first, and rather uncontrollable, approximation for the re¬ 
alistic tight-binding models of Weyl semimetals. First of 
all, the introduction of a momentum cutoff for a contin¬ 
uum Dirac Hamiltonian violates the vector current con¬ 
servation, which should then be manually restored by 
adding the so-called Bardeen counterterm [18, 27]. Sec¬ 
ond, the continuum Dirac Hamiltonian possesses exact 
chiral symmetry, which is typically violated in lattice sys¬ 
tems at energies sufficiently higher than the energies of 
the Dirac points. This violation drastically changes the 
pattern of spontaneous symmetry breaking, in particu¬ 
lar, turning the phase with spontaneously broken con¬ 
tinuous U (1) chiral symmetry and massless Goldstone 
modes into the “axionic insulator” phase (or the Aoki 
phase in the terminology of lattice QGD) in which only 
discrete Z 2 CV symmetry is spontaneously broken and no 
Goldstones emerge [41, 42]. Therefore it is important to 
go beyond the Dirac cone approximation and to perform 
a consistent calculation of the chiral magnetic conductiv¬ 
ity directly in the realistic tight-binding model of a Weyl 
semimetal. 

In this paper we follow [17, 43-49] and use the Wilson- 
Dirac Hamiltonian as the simplest realistic lattice model 
of Dirac semimetals and/or topological insulators, which 
are then tuned into a Weyl semimetal phase by adding 
the parity- or time-reversal-breaking terms. In our case, 
parity is broken by the chiral chemical potential ^j^a- We 
will use the mean-field approximation both to find the 
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ground state of our model and to calculate the static 
chiral magnetic conductivity. Since the static chiral mag¬ 
netic conductivity is the static linear response of the elec¬ 
tric current to a small static external magnetic field, in 
the mean-field approximation it can be related to the lin¬ 
earized response of the static mean-field condensates to 
the static external gauge field. Having in mind that re¬ 
cent experiments [4—6] are performed at very low temper¬ 
atures (~ 10^ K), we also mostly concentrate on the case 
of zero temperature, and discuss the effect of finite tem¬ 
perature (which in most cases acts against CME) only 
briefly. 

In contrast to the results obtained in [27] for the con¬ 
tinuum Dirac Hamiltonian, we find that the only source 
of enhancement of the chiral magnetic conductivity is the 
renormalization of the chiral chemical potential and 
it never exceeds the universal value acME = ^ 2 -^^ with 
the renormalized fj,A- The universal value (2) is only ap¬ 
proached in the regions of the phase diagram where our 
model features massless Dirac cones in the energy spec¬ 
trum. We find that in these regions the corrections due to 
inter-electron interactions (apart from those originating 
from the renormalization of /ta) are very small, not larger 
than a few tenth of percent. Corrections due to interac¬ 
tions only become important when the energy spectrum 
has a large gap, and reach maximum in the phase with 
spontaneously broken CV symmetry. 

On the other hand, the phase structure of our lattice 
model in the parameter space of the chiral chemical po¬ 
tential and interaction strength is more similar to the one 
obtained for the continuum Dirac Hamiltonian [27]. In 
particular, we find that interactions effectively enhance 
the chiral chemical potential ^a- Nonzero chirality im¬ 
balance also shifts the boundaries of the Aoki phase to¬ 
wards weaker inter-electron interactions - much like the 
critical interaction strength for spontaneous breaking of 
chiral symmetry is decreased in the presence of ^a [27]. 
The only qualitative difference with the continuum case is 
that in the lattice model we do not observe any Cooper- 
type instability towards spontaneous breaking of chiral 
symmetry at small inter-electron interactions. The ori¬ 
gin of this discrepancy is quite clear: for the Wilson-Dirac 
Hamiltonian, continuous chiral symmetry is already ex¬ 
plicitly broken at high energies. Preliminary mean-field 
studies of the phase diagram of Wilson-Dirac fermions 
with chirality imbalance and on-site inter-electron inter¬ 
actions were already reported by two of the authors in 
[42]. 

The outline of this paper is the following: in Section 2, 
we present a general framework for the mean-field calcu¬ 
lations of the phase diagram and of the chiral magnetic 
conductivity in our lattice model. In Section 3, we study 
the mean-field phase diagram of the Wilson-Dirac Hamil¬ 
tonian with chirality imbalance in the parameter space 
of chiral chemical potential, bare Dirac mass and inter¬ 
electron interaction potential. In Section 4, we perform 
the mean-field calculation of the static chiral magnetic 
conductivity in different regions of the phase diagram. 


In Section 5 we comment on the feasibility of the ex¬ 
perimental measurements of the static chiral magnetic 
conductivity. We summarize the results of our study in 
the concluding Section 6 . 

2. MEAN-FIELD APPROXIMATION FOR THE 
WILSON-DIRAC HAMILTONIAN WITH 
ON-SITE INTER-ELECTRON INTERACTIONS 

The starting point of our study is the Wilson-Dirac 
Hamiltonian with chiral chemical potential and on¬ 
site inter-electron interactions with repulsive potential 
C7 > 0: 

= E E - 2) (4) 

x,y X 

where V’J and are the fermionic creation and annihi¬ 
lation operators at the lattice sites x and y and h^xy is 
the single-particle Wilson-Dirac Hamiltonian on the cu¬ 
bic three-dimensional lattice. In the absence of external 
fields, it has the following form in momentum space: 

3 

= E aiVF sin (ki) + 

3 

+2r 7o E + 70™^°^ -k 75 /i^^, (5) 

i=l 

where vr is the Fermi velocity, r is the Wilson coefficient, 
is the Dirac mass term, 

Ofc = diag (cTfe, -cTfe), 75 = diag (I, -I ), 



are the Dirac matrices, crfc, k = I, 2, 3 are the Pauli ma¬ 
trices and / is the 2x2 identity matrix. In (5), we have 
expressed all dimensionful quantities in units of the in¬ 
verse lattice spacing a~^, which is hence just unity in 
what follows. We also use the units with h — c = 1. 
The single-particle Hamiltonian (5) is the simplest model 
which mimics most important features of the more real¬ 
istic tight-binding models of Dirac and Weyl semimet¬ 
als [17, 43-49], such as the emergent chiral symmetry at 
small energies. 

The explicit coordinate-space form of y in (4) with 
external gauge fields A {x) can be found by the Peierls 

T + fc 

substitution of the link phases Ax^k = / dzAk ( 2 ) into 

X 

the finite-difference operators on the lattice: 

3 

^x,y ~ ^ ^ ictfc'CpV k xy ^ ^xy 

k^l 

-k 7 om(°) +75 Ma\ 


(7) 
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where 

Vfcxi/ = 2 (® '°’'°^x+k,y ~ ^ ^ 

A^y = E (2 - - ^-^^-^■^Sx-k,y) (8) 

k=l 

are the lattice discretizations of the covariant derivative 
operator and (minus) the Laplacian operator — 
and fc is a unit lattice vector in the direction k. 

In (4) we include only instantaneous on-site inter¬ 
electron interactions. On the one hand, mean-field phase 
diagrams obtained with only on-site interactions in most 
cases qualitatively reproduce most important features of 
the phase diagrams with more realistic interaction po¬ 
tentials. On the other hand, the Hubbard-Stratonovich 
transformation which is an important ingredient of the 
mean-field analysis greatly simplifies for on-site interac¬ 
tions. We also neglect the magnetic part of electromag¬ 
netic interactions, which is suppressed as 1/vp, and treat 
the vector gauge potential {x) in the Hamiltonian (4) 
as a non-dynamical external field which only serves as a 
source for the static electric current. The instantaneous 
approximation for our interaction potential can also be 
justified by the smallness of the Fermi velocity vp in re¬ 
alistic systems. 

Moreover, as we show in Appendix B, if we only con¬ 
sider instantaneous inter-fermion interactions, the depen¬ 
dence on the Fermi velocity amounts to simple linear scal¬ 
ing of the model parameters and observables. Therefore 
we set vp = 1 and summarize the most important formu¬ 
las from which the full dependence on uf can be restored 
in all our results in Appendix B. The dependence on the 
Wilson coefficient r in (7), (5) is strictly speaking math¬ 
ematically nontrivial. However, in practice changing r 
results only in a rescaling of the phase diagram in the 
direction of the bare Dirac mass [41], and no qual¬ 
itative changes of the phase structure occur. For this 
reason we also set r = vp = 1 in order to reduce the 
number of free parameters in our model and simplify the 
analysis. 

Below we summarize the basic steps leading to the 
mean-field approximation for the partition function with 
the Hamiltonian (4). A more detailed derivation can be 
found in [27]. In order to calculate the partition function 

2^ = Tr exp we perform the standard Suzuki- 

Trotter decomposition of the exponent exp(^—H/T^, 

splitting the Euclidean time r G [O ... T“^] into infinites¬ 
imal intervals of size At: 


use the Hubbard-Stratonovich transformation to rewrite 
the exponents involving the interaction term as 

exp (^-At Hi"j = 

exp ^-At^ ( '^^4^*''^^ ( 10 ) 

where (t) is a Hermitian traceless matrix with two 
Dirac spinor indices [63] and we include some trivial nor¬ 
malization factors into the definition of the integration 
measure d^x {t). In the following it will also be con¬ 
venient to represent $ 3 , as a sum over a basis set of 15 
traceless Hermitian spinor matrices: 

15 

^x = Y^ ^A<^x,A, 

A^l 

Ta = {75Q;fc,7o,7o75afe,-*757o,*7oafc,75,a/c}, (H) 

where k = 1,2,3 labels spatial directions. The ma¬ 
trices Fa are normalized as Tr (FaFs) = 4i5a_b, so 
that the action of the Hubbard-Stratonovich field reads 

Tr^l/ m = ^IJU. 

After inserting the transformation (10) into the 
Suzuki-Trotter decomposition (9), we can represent the 
partition function Z = Tr exp in terms of the 

partition function Z [<I> 2 , (r)] of a free fermion gas in 
the background of the space- and time-dependent field 
(r): 


Z = y (r) Z [$, (t)] X 


X exp — / dr 


Tr4>g (r) 

417 

(r)] = 


= TrT exp 


- I + ’ (12) 


with the time ordering operator T. 

We now apply the mean-field approximation and re¬ 
place the integral over (t) by its value <i)J (r) at the 
saddle point of the path integral (12). Moreover, we as¬ 
sume that $* (r) = does not depend on the Euclidean 
time T. The free energy T = —TlnZ which corresponds 
to the partition function ( 12 ) is then simply the func¬ 
tional 


Tr exp = 

= ^lim^ Tr ^ 

where Hq is the free part of the Hamiltonian (4) and 
•il^x — 2^ is the interaction term. Next we 


^[<i>.] = J-o[d>.] + 5:^, (13) 

X 

where one should take the value $ 3 , = $* which min¬ 
imizes J- [ 4 ) 3 ;] and J-Q [ 4 ) 3 ;] is the free energy of a free 
fermion gas with an effective single-particle Hamiltonian 

hxy = h^xy + ^xdxy (14) 
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The explicit expression for Fq reads 

-Fo [$x] = -T ^ In (l + , (15) 


where are the energy levels of the single-particle Hamil¬ 
tonian (14). 

The static electric current which enters the definition 
(2) of the chiral magnetic conductivity can be expressed 
as the variation of the free energy F over the gauge po¬ 
tential Ak (x). More precisely, in the lattice model (4) 
the operator of electric current density is associated with 
lattice links, and its expectation value is related to the 
derivative of the partition function over the link phase 
factors Ax^k- {jx,k) = Obviously, this current 

should vanish in the absence of external gauge fields. The 
linearized response of the electric current to a small ex¬ 
ternal gauge field is then given by the second variation of 
the free energy with respect to the vector gauge potential 
at zero external field: 


{jx,k ) 'y ' 

v,l 


S^J^[A,,k] 

^Ax^k ^Ay l 


A. 




A=o 


(16) 


where we have explicitly restored the dependence of the 
free energy on the link phases Ax^k- After expressing the 
static external magnetic field as Bi (x) = Cijk dxj (^) 
and rewriting the relation (16) in the momentum space, 
we arrive at the following Kubo relation for the static 
chiral magnetic conductivity (2) as a function of the wave 
vector k of an external magnetic field [50-52]; 

crcME ih) = = 

^ x,V 

=_L V Afcafe-y.d S'^J^[Ax,k] 

*3 ^Ax^i 6Ay^2 


(17) 


For definiteness, we have assumed in the above expression 
that the wave vector of the external magnetic field is 
parallel to the third coordinate axis and the vector gauge 
potential is parallel to the second axis. In this case the 
magnetic field and the CME current are parallel to the 
first coordinate axis. 

In the mean-field approximation, we replace the full 
free energy F = —TlnZ by the minimal value of the 
mean-field free energy (13). When using (17) to calcu¬ 
late the chiral magnetic conductivity, one should now 
keep in mind that the extremum value of the Hubbard 
field = $* \Axa] which minimizes the free energy 
F[^xj-^x,k\ given by (13) in general depends on the ex¬ 
ternal gauge field Ax k ■ The variation of the 

extremum value [Ax^k] of the Hubbard-Stratonovich 
field over the external gauge field Ax^k can be found as 
the derivative of an implicit function using the extremum 
equation 


dF[^x.AxA 


9$ 


x,A 


= 0 . 


(18) 


After some algebraic manipulations which are summa¬ 
rized in detail in [27], we find the following general ex¬ 
pression for the second derivative over the vector gauge 
potential which enters (17): 


5^F[^l,AxA d^F[^l,AxA 


- E G 


z.A’.t.B 


d^F[^x,Ax,k] 

dAx^r 


dAyj d^t,B 


,(19) 


where Gz,A-,t,B is the propagator of the Hubbard- 
Stratonovich field defined as: 


E/ 

y,B 


d^F[^x,Ax,k] 

d^y^Bd^Z,C 


^xzSaC- 


( 20 ) 


Now the derivatives djdAx^i and d/d^x,A in (19) and 
(20) should be calculated as partial derivatives, that is, 
in the process of differentiation the mean-field free energy 
functional (13) should be considered as a functional of the 
two independent fields Ax^i and ^x,A- 

The first term in (19) describes the electromagnetic re¬ 
sponse of the non-interacting system with renormalized 
single-particle Hamiltonian (14). Since the parameters 
of the renormalized, rather than the bare, single-particle 
Hamiltonian are physically observable, this term does not 
really describe the effect of interactions. Only the sec¬ 
ond term in (19) describes nontrivial corrections to the 
current-current correlators due to inter-electron interac¬ 
tions. In the weak-coupling regime, it can be represented 
as an infinite sum of ladder diagrams [27]. 

While the differentiation of the second term in (13) 
is trivial, the calculation of the second derivatives 






dA^^idAyj’’ dA^_idiy,A 


and 




of the fermionic 


free energy (15) is somewhat lengthy, and we summa¬ 
rize it in Appendix A. The basic idea of this calculation 
is that in the limit of static Hubbard-Stratonovich and 
vector gauge fields the derivatives of the single-particle 
eigen-energies which enter (15) can be obtained simply 
using the second-order quantum mechanical perturbation 
theory for the effective single-particle Hamiltonian (14): 


= (*.| |».) + 


80 d£, 


8^89 




('I'd t i I'f'.) 


E 


ei - e, 


( 21 ) 


where the variables 9, ^ can be either the link factors 
Ax,i or the Hubbard-Stratonovich fields ^x,A and |'I'i) 
denotes the eigenstate of the single-particle Hamiltonian 
(7) which corresponds to the energy level e^. Note also 
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that since the link factors enter the single-particle 
Hamiltonian (7) in a nonlinear way, we should in gen- 
eral keep the terms involving its second derivative 
These terms are absent for the continuum Dirac Hamil¬ 
tonian [27], where the gauge vector potential enters lin¬ 
early, essentially in the same way as the component 
of the Hubbard-Stratonovich field which corresponds to 
Fji = Oi- It is precisely this difference which renders 
the current-current correlator in (16) manifestly gauge- 
invariant and effectively replaces the subtraction of the 
Bardeen counterterm. 

Differentiating now the free energy of the free fermion 
gas (15) and using the above expressions, we obtain 


92 


ddd£ ^ 

i 

= E 


9^/(ei) da dti ^ df{ei) 
9e2 de d^ 

-1 


de do d^ 


i ^ 

fl^.) 


Ci - e, 


-^n(e,) 
i¥=j 

, V- . Ill'll*) 

+ > n (Ci)----,(22) 

e* — £7 


where we have denoted /(e) = —Tin (l-|-and 
n (e) = df (e) jde = is the Fermi factor. 

Since in (17) we are calculating acME in the linear 
response regime ( 22 ) in the absence of background mag¬ 
netic fields, in the expression ( 22 ) we should set to 
zero and the Hubbard-Stratonovich field $ 2 , - to its mean- 
field value. In this paper, we assume that the mean-field 
configuration of the Hubbard-Stratonovich field is spa¬ 
tially homogeneous (in the following we discuss the va¬ 
lidity of this assumption in more details). In this case, 
the eigenstates j'l'i) can be found explicitly, which re¬ 
duces the calculation of the second derivatives ( 22 ) which 
enter (19) to a simple numerical summation over lattice 
momenta. In Section 3 we present the mean-field phase 
diagram of our model (4). In Section 4 we then use the 
mean-field single-particle Hamiltonian to calculate the 
chiral magnetic conductivity (17) using the expressions 
(19) and (22). 


3. MEAN-FIELD PHASE DIAGRAM IN THE 
PRESENCE OF CHIRALITY IMBALANCE 

To compute the mean-field phase diagram of the 
model (4) we assume that the saddle-point Hubbard- 
Stratonovich field = 4>* is homogeneous and does 
not break rotational symmetry. Under this assumption 
the saddle-point values are restricted to the form 

$* = (mr - 7o -I- imao-f5 + 75,(23) 


where we introduced the renormalized mass rrir, the 
renormalized chiral chemical potential fj-A and the CV- 
breaking mass term rrii which corresponds to the “axion” 
condensate (' 0 ^ 7 o 75 ' 0 )- With the mean-field value (23) 
for the Hubbard-Stratonovich field and in the limit of 
vanishing temperature the free energy density becomes 


T[^*] 


1 

if 


£<0 


(777(0) _ ^^^2 _|_ 7772 _|_ _ ^^^2 


(24) 


The first term is a sum over all negative energy levels 


e = £s,a = s\J{S - a^iA^ + rnl + (m^ -I- Wf (25) 

of the effective single particle Hamiltonian (14) and for 
the sake of a compact notation we define s = ± 1 ,ct = ± 1 , 
S = sin^(ki) and W = ^ 2r sin^(fci/2). 

V ki ki 

Before discussing the phase diagram of our model, it 
is useful to briefly recollect some basic facts about the 
Wilson-Dirac Hamiltonian (5). In the context of lattice 
QCD Wilson added the term proportional to A^y to get 
rid of the so-called “doublers”, unwanted fermionic de¬ 
grees of freedom that should not be present in the contin¬ 
uum limit of lattice QCD. By adding the Wilson term to 
the naive discretization the doublers acquire a mass that 
depends on the inverse lattice spacing. In the continuum 
limit the doublers become infinitely heavy and decouple 
from the theory. We use Wilson-Dirac fermions to model 
the band structure in a crystal and use a fixed lattice 
spacing that was set to one in equation ( 8 ). In this con¬ 
text, the Wilson term A^y reflects the fact that in Dirac 
or Weyl semimetals the continuum chiral symmetry is 
only an emergent symmetry at low energies, which is ex¬ 
plicitly broken away from the Dirac/Weyl points. The 
mass of the doubler fermions is given by the poles of the 
propagator in momentum space, which are located at po¬ 
sitions where ^ 7 ,. sin(fc 7 )^ = 0. Since the momenta lie 
in the Brillouin zone, i.e. ki S (—tt, tt], the mass of the 
doublers reads 


rrid = nir + 2A, (26) 

where A is the number of momentum components ki 
with ki = TT and we assume that the Wilson parameter 
r = 1 according to the discussion in the previous Section. 
In three dimensions A can take the values A = 0,1,2,3. 
There is only one way to achieve A = 0 and A = 3, all 
entries of k have to be the same. For A = 1 we have three 
possibilities to place the entry tt in a three dimensional 
vector, 

fc = (0, 0, tt)"^ , fc = (0, TT, 0)'^ and k = (tt, 0, O)^ . (27) 

Analogous for A = 2 there are three ways to place the 0. 
This means that for a bare mass of nir = 0 and = — 6 
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we have Nf = 1 Dirac cones in the energy spectrum of 
our model, whereas there are Nf = 3 Dirac cones for 
rrir = —2 and = —4. It follows from equation (26) 
that non-trivial phases with massless excitations only ex¬ 
ist for —6 < 771^°^ < 0 and we restrict our study to this 
parameter range. 

After setting r = uf = 1 and expressing all dimen¬ 
sionful quantities in units of the lattice spacing we are 
thus left with three free parameters in our model: the 
bare mass the inter-fermion interaction potential 

U and the bare chiral chemical potential To com¬ 
pute the phase diagram we scan over these parameters 
and find the renormalized values of the parameters 
rrii and of our model by numerically minimizing the 
free energy (24) with respect to them. 



FIG. 1: The CP-breaking mass term rrii in (23) as a function 
of interaction potential U for different values of the bare chiral 
chemical potential 

Our model exhibits two distinct phases: A CV- 
symmetric phase with rrii = 0 and the “Aoki phase” 
where the CV symmetry is spontaneously broken by a 
nonzero CT^-breaking effective mass term rrii ^ 0. In our 
lattice model, this spontaneous breaking of a discrete CV 
symmetry is a remainder of the spontaneous breaking of 
continuous chiral symmetry which is now only an emer¬ 
gent low-energy symmetry. In the context of condensed 
matter physics, the Aoki phase is the “axionic insulator” 
phase with a nonzero condensate of an effective axion 
field (see e.g. [53-55] for some examples of effective ax¬ 
ion fields in condensed matter systems). On Fig. 1 we 
show the mean-field value of the CT^-breaking mass term 
rrii as a function of interaction potential U for different 
values of the bare chiral chemical potential /r®. For all 
values of the chiral chemical potential the transition 
between the normal insulator/semimetal phase and the 
Aoki phase with ^ 0 appears to be a sharp second- 
order phase transition. 

We plot the mean-held phase diagram of our model 
in the m*^°^-[/-plane for different values of the bare chi¬ 
ral chemical potential on Fig. 2. The phase diagram 
is symmetric around the line of = —3.00. Increas- 



FIG. 2: Mean-field phase diagram of our model in the 
plane for different values of the bare chiral chemical potential 
fj,A- The points mark the boundaries of the Aoki phase. 

ing the bare chiral chemical potential /x^^ shifts the on¬ 
set of the Aoki phase to slightly smaller values of the 
inter-electron interaction U. A similar effect of the chiral 
chemical potential on the critical interaction potential for 
spontaneous breaking of chiral symmetry has been found 
in [27] for continuum Dirac fermions. 

In a previous study [42] we presented results for 
Ls = 10, since numerical calculations for Lg = 8,10,12 
indicated that finite size effects have little influence on 
the mean-field phase diagram. This holds true in general 
but we observe that the thin structures which point to 
the characteristic values of = 0, —2, —4, —6 and are 
often called “Aoki fingers” get thinner as Lg is increased. 
In the case of Lg = 50 we were no longer able to resolve 
the fingers numerically, see Fig. 13. Even for small lat¬ 
tice sizes the fingers are only present at zero bare chiral 
chemical potential and vanish as soon as ^ 0. For 
the calculation of the chiral magnetic conductivity we are 
interested in values 0 < < 0.3, where the phase dia¬ 

grams for Lg = 10 and Lg = 50 lie on top of each other. 
Therefore the volume dependence of the Aoki fingers at 
= 0 is of little interest for the main points of this 
paper and we discuss the origin of finite size effects in 
Appendix C. 

Similarly to the case of continuum Dirac fermions [27], 
our lattice model also exhibits a strong renormalization 
of the chiral chemical potential. Figure 3 shows the chi¬ 
ral chemical potential as a function of the inter-electron 
interaction U for different values of the bare chiral chem¬ 
ical potential /x^^ We note that the growth of fj,A with 
U becomes somewhat more pronounced within the phase 
with spontaneously broken CV symmetry. 

Let us now discuss in more detail the validity of our 
assumption of the spatial homogeneity of the fermionic 
condensates. In the process of calculating the current- 
current correlator (19) we have to compute the Hes¬ 
sian matrix of second derivatives of the 

mean-held free energy with respect to the Hubbard- 
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FIG. 3: Renormalized chiral chemical potential ha as a func¬ 
tion of the interaction potential U for different values of the 
bare chiral chemical potential ■ 
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FIG. 4: Phase diagram in the parameter space of bare mass 
and the interaction potential U at different tempera¬ 
tures. The chiral chemical potential is Hxed at =0.1. 


Stratonovich field ^x,A in the vicinity of the homoge¬ 
neous configuration $2; = $* which minimizes the mean- 
field free energy functional (13). If the homogeneous con¬ 
figuration $2^ = $* is a local minimum in the space of 
Hubbard-Stratonovich fields ^x,A, all the eigenvalues of 
this matrix should be positive. Since $* is a spatially 
homogeneous field configuration, the entries of this Hes¬ 
sian matrix depend only on x — y, and we can calcu¬ 
late its eigenvalues by performing the Fourier transfor¬ 
mation with respect to x — y. At each wave vector in 
the Fourier transform, one then has to diagonalize only 
a 15 X 15 matrix with entries corresponding to differ¬ 
ent spinor structures of $2: (see (11)). We have calcu¬ 
lated the eigenvalues of the Hessian matrix of T [<&] for 
every point in the parameter space which we have con¬ 
sidered in our mean-field calculations and for wave vec¬ 
tors fca = 2TTm/Ls, m S {0,..., 8} and found that they 
are always positive. This calculation suggests that the 
uniform saddle-point configuration $* of the Hubbard- 
Stratonovich field defined in (23) is indeed a local min¬ 
imum in the space of all possible field configurations of 
<i>2; and if there is a non-homogeneous field configuration 
which further decreases the free energy it should be sep¬ 
arated from the homogeneous configuration $* by some 
potential barrier. These findings to some extend justify 
the assumption of homogeneity made for the saddle point 
value of the field. In the concluding Section 6, we fur¬ 
ther discuss the plausibility of this assumption in view of 
some recent theoretical studies which predict inhomoge¬ 
neous ground state for chirally imbalanced fermions. 

Finally, we briefly consider the effect of finite temper¬ 
ature on the phase structure of our model. On Fig. 4 
we plot the boundaries of the Aoki phase in the parame¬ 
ter space of bare mass and the interaction potential 
U at different temperatures. The bare chiral chemical 
potential is fixed at =0.1. As expected, finite tem¬ 
perature tends to destroy the fermionic condensates and 
hence shifts the phase with spontaneously broken CV 


symmetry towards stronger interactions. This shift is, 
however, quite small even for temperature T = 1 com¬ 
parable to the band width in our model. We conclude 
therefore that the effect of temperature is rather small at 
the level of the phase diagram. 


4. MEAN-FIELD CALCULATION OF THE 
STATIC CHIRAL MAGNETIC CONDUCTIVITY 

Once we have found the renormalized values of nir, rrii 
and Ha for our model, they can be used to calculate the 
chiral magnetic conductivity <JcME{k) from equations 
(19) and (17). On Fig. 5 we plot ctcme (k) as a function 
of the wave vector k (which is parallel to one of the co¬ 
ordinate axes) of the external magnetic field at several 
characteristic points in the phase diagram of our model. 
The bare chiral chemical potential is Ha'^ = 0-05 every¬ 
where. Within each plot, we scan the phase diagram 
along the line of constant renormalized mass mr- Since 
for Dirac/Weyl semimetals one is typically interested in 
the regions with gapless spectrum, we have considered 
the values uir = 0, nir = —2 and nir = —4. According 
to the discussion in Section 3 above, at = —2, —4 
and rrir = 0, —6 the dispersion relation of our model has 
Nf = 3 or Nf = l Dirac cones, correspondingly. This 
results in an obvious linear scaling of acME (k) with Ay, 
which can be clearly seen on Fig. 5. Moreover, one can 
show that the renormalized masses rrir and —6 — rrir are 
equivalent, but the roles of left- and right-handed spinor 
components interchange between them. In other words, 
the chiral chemical potential effectively changes sign be¬ 
tween rrir and —6 — m^. This can be clearly seen on 
Fig. 5 for the lines with rrir = ~2 and rrir = ~4, for 
which acME (k) only differs by a sign. We note also that 
by virtue of such change of sign between the values rrir 
and —6 — nir the chiral magnetic conductivity acME {k) 
exactly vanishes for all k on the line rrir = ~3, which is 
the symmetry axis of the phase diagram of our model. 
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FIG. 5: Numerical results for the chiral magnetic conductivity ucme (k) at characteristic values of nir for different values of 
the interaction potential U (lines with points). Clockwise from top left: = 0.00, = —2.00,mr. = —4.00 and rrir = —2.20. 

In all cases = 0.05 and Ls = 50. For comparison we also plot the result (28) for free continuum Dirac fermions with 
Pauli-Villars regularization and with the renormalized chiral chemical potential fiA (lines of the same style as for numerical 
data, but with no points on them). 


With these model-dependent features taken into ac¬ 
count, we observe a qualitative agreement between our 
data at small k outside of the Aoki phase and the con¬ 
tinuum result for (Jcme (k) obtained with Pauli-Villars 
regularization [ 20 , 21 ]: 


{2Ef 


k'A ~ 1 

MA + ^ , log 


^CME (^) — 

— k 


2^a “t" k 


(28) 


In particular, one can see the characteristic growth of 
crcME (k) around k = 2fiA- At fc > 2fiA ctcme {k) ap¬ 
proaches the asymptotic value ± 5 ^ of the continuum 
result, but at even larger k we observe a decrease due to 
explicit chiral symmetry breaking at larger momenta. It 
is important to stress that the behavior of ocme (k) de¬ 
pends on the renormalized chiral chemical potential ^A^ 
which is an experimentally measurable quantity. E.g. the 
angular-resolved photoemission spectroscopy (ARPES) 
measurements in the vicinity of the Dirac cones (see for 
example [5]) should yield the renormalized chiral chemi¬ 
cal potential ^a- Therefore we use the renormalized value 


of fiA in order to compare our data with the continuum 
result (28). 

We see that for all plots the lattice result for acME {k) 
is smaller than the continuum result (28) calculated with 
the renormalized chiral chemical potential. We conclude 
therefore that the observed enhancement of ucme (k) 
is mostly due to the growth of fj-A with interaction 
strength. As the effective single-particle Hamiltonian be¬ 
comes gapped due to a nonzero C'P-breaking mass term 
in the Aoki phase, acME {k) becomes significantly sup¬ 
pressed at all k, despite the even stronger growth of ^a 
in this phase. In order to illustrate the behaviour of 
crcME (k) in the case when the spectrum of the effec¬ 
tive single-particle Hamiltonian is gapped outside of Aoki 
phase, on Fig. 5 (at the bottom on the left) we also plot 
it along the line with = —2.20. We again see that the 
CME response becomes suppressed, as could be expected 
also from continuum calculations [ 20 , 21 ]. 

Let us now check whether the behavior which we have 
found for the above specific values of model parame¬ 
ters is typical throughout the whole phase diagram of 
our model. To this end, it is advantageous to char- 
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acterize (Jcme (k) in terms of its asymptotic value at 
/iA ^ fc ^ 1. In this range of momenta the lattice re¬ 
sult for (JcmeO^) should approach the continuum one, 
which tends to the asymptotic value at k 

(see Fig. 5). The universality of this asymptotic value is 
due to its relation to the anomaly coefficient, which is, 
however, quite nontrivial [20]. 

However, our calculations are performed on a finite 
lattice, and the values of k are quantized. Therefore we 
have a rather narrow range of lattice momenta fiA 
fc 1 for which it makes sense to compare the chiral 
magnetic conductivity crcME{k) with have 

chosen to estimate the asymptotic value a (k) at ^a ^ 
fc <C 1 simply as the maximal value of acME (k) among 
all discrete momenta which we consider. We denote this 
estimate as cTmax'- 

^max = max acME (k) ■ (29) 



FIG. 7: Numerical estimate Umax of the asymptotic values of 
ucME (fc) at /iA <C fc <C 1 as a function of the renormalized 
chiral chemical potential /ta for different values of the bare 
mass and the interaction potential. The points are connected 
with straight lines to guide the eye. For small pA Umax is a 
linear function of the renormalized chiral chemical potential. 
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FIG. 6: Numerical estimate Umax of the asymptotic value of 
ucME (fc) at PA <C fc <C 1 as a function of bare mass 
and the interaction potential U at fixed bare chiral chemical 
potential p® = 0.05. The red crosses mark the border of the 
Aoki phase and the black points mark the lines of constant 
nir = 0, —2, —4, —6. 

On Fig. 6 we plot our numerical estimate Umax of 
o'CME (k) at PA ^ ^ I in the parameter space of 

the bare mass and the interaction potential U for 
~ 0-05. In agreement with our previous observa¬ 
tions, we find that Umax is peaked around the lines of con¬ 
stant renormalized mass rrij. = 0, —2, —4, —6, at which 
our model has massless excitations. In the gapped phase 
with spontaneously broken CV symmetry, Umax quickly 
decreases. In general, the behavior of Umax is non-trivial 
and depending on the location on the bare mass axis 
increasing the inter-electron interaction strength U can 
lead to an increase or a decrease of Umax- 

In order to compare our estimate of the asymptotic 
values of ucme (k) at pA ^ fc ^ I with the universal 
value in (2) which is related to axial anomaly coef¬ 
ficient, we further assume that Umax is a linear function 
of PA- In order to justify this assumption, on Fig. 7 we 
show the dependence of Umax on pA for several different 
points in the parameter space of our model. We thus 


estimate the ratio of the asymptotic value of ucme (k) 
at Pa ^ fc ^ I to the renormalized chiral chemical po¬ 
tential PA by fitting Umax at p^^ = 0.05, O.IO, 0.15 with 
a linear function Umax = A^la, where A is the fitting 
parameter. We note also that in realistic experiments, it 
is difficult to achieve large values of pA comparable to 
the band width, so in practice one is anyway interested 
in small values of pA- 



FIG. 8: The ratio of the estimated asymptotic value Umax of 
UCME (fc) at fiA <C fc <C 1 to the renormalized chiral chemi¬ 
cal potential fiA- The black dots mark the lines of constant 
renormalized mass = —6, —4, —2 and nir = 0, respec¬ 
tively, and the red crosses indicate the border of the Aoki 
phase for p)]"* = 0.0. 


Our numerical estimate of the ratio at 

fJ-A MA 

PA ^ ^ I is shown on Fig. 8 as a function of the 

bare mass and the interaction potential U. Just as 
on Fig. 6 we see that this ratio is peaked along the lines 
of constant rur = 0, —2, —4 and —6, where the fermionic 
single-particle spectrum is gapless, and is strongly sup¬ 
pressed in the vicinity of the Aoki phase. We see again 
that even in the vicinity of the lines with Nf = I and 
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Nf = 3 massless Dirac cones in the energy spectrum the 
ratio never exceeds the universal value TV// (27r^). 
We thus conclude that strong enhancement of acM e ob¬ 
served for the continuum Dirac fermions with sponta¬ 
neously broken chiral symmetry in [27] does not happen 
in the lattice model with only emergent chiral symmetry. 

It is also interesting to check what is the nontrivial 
contribution of inter-electron interactions to <7cme {k). 
As discussed in Section 2, in the general expression (19) 
for the current-current correlators the first, tree-level, 
summand describes the change of the electromagnetic re¬ 
sponse merely due to renormalization of the parameters 
rtir and of the effective single-particle Hamilto¬ 
nian, whereas the second summand describes nontrivial 
loop corrections originating from ladder diagrams with 
arbitrary number of fermionic loops. In order to quan¬ 
tify the importance of interactions, on Fig. 9 we sepa¬ 
rately plot the tree diagram contribution to acME {k) 
(first summand in (19)) and the contribution of loop di¬ 
agrams (second summand in (19)) as a function of fca for 
TOr = 0, = 0.05 and vrir = —2, = 0.05. We see 

that the loop contributions are numerically very small in 
the phase with unbroken CV symmetry - at maximum 
around one percent for rrij. = —2, and by almost two or¬ 
ders of magnitude less for rur = 0. It is interesting to 
note here the dependence on the number of Dirac cones 
in the spectrum. The only exception is the point with 
rUr = — 2 and U = 2.25 in the phase with broken CV 
symmetry, where the contribution of the tree diagram 
becomes very small due to gap opening, but the loop 
contributions stay almost the same as in the unbroken 
phase, thus they become more important. 

We further illustrate the importance of loop contribu¬ 
tions on Fig. 10 by plotting the ratio of the tree diagram 
and the loop contributions to (Jcme (k) at fixed k = 6 
as a function of the bare mass and the interaction 
potential U. We see again that the loop contribution is 
very small along the lines with massless Dirac cones in 
the spectrum, and only becomes important in the regions 
of phase diagram where the energy spectrum has a large 
gap and thus the tree diagram contribution is strongly 
suppressed. Within the “axionic insulator” phase with 
broken CV symmetry, the importance of loop contribu¬ 
tions reaches its maximum. In particular, it is interest¬ 
ing to note that the loop contribution is strongly peaked 
along the two lines starting roughly from = —5 and 
= — 1 and going into the region of the broken CV 
phase. 

Finally, on Fig. 11 we illustrate the effect of finite tem¬ 
perature on the chiral magnetic conductivity (Jcme (k). 
For all values of the model parameters finite tempera¬ 
ture suppresses the CME response. The suppression is 
not very large if the temperature is comparable with the 
chiral chemical potential. At higher temperatures of the 
order of the band width, acME (k) becomes almost zero. 
Such temperature dependence of the chiral magnetic con¬ 
ductivity is expectable, since the temperature tends to 
reduce the differences in the occupation numbers of left- 


and right-handed fermions of the same momentum. In¬ 
deed, in the limit of infinite temperature, the Fermi dis¬ 
tributions in (15) and (22) become independent of /i^, 
and hence the CME response should vanish. 

5. CONNECTION WITH EXPERIMENT 

Spatially modulated magnetic fields are very conve¬ 
nient for theoretical considerations of the chiral magnetic 
effect, in particular, for establishing its connection with 
the anomaly coefficient. In experiments, however, one 
typically subjects the Weyl semimetal sample to a static 
and constant magnetic held. Since our calculations show 
that cfcME{k) vanishes in the limit fc —?► 0, a natural 
question is how to measure such a quantity using a con¬ 
stant magnetic held. While the static Kubo formula (17) 
on which we base our calculations probably does not take 
into account all the details of the realistic experimental 
setup, in this Section we show that nevertheless our re¬ 
sults provide a reasonable estimate of the CME response 
in experiments, at least by order of magnitude. 

To this end we note that the real samples of 
Weyl/Dirac semimetals used in experiments have hnite 
size of order of 1mm [5, 6]. Correspondingly, one can ex¬ 
pect that magnetic helds have spatial modulation with 
characteristic wavelength of the same order or even 
shorter, due to inhomogeneities and grain structure of 
real crystals. 

In this paper, we have performed calculations for hnite- 
size systems with periodic boundary conditions, which 
do not allow for arbitrarily small constant magnetic 
fields due to magnetic flux quantization constraints. On 
the other hand, outside of the axionic insulator phase 
and for sufficiently small momenta our results are quite 
well described by the expression (28) obtained for Dirac 
fermions with Pauli-Villars regularization in infinite con¬ 
tinuum space. In order to model the spatial modulation 
of magnetic field due to the finite sample size in a way 
which is consistent with the approximations used in this 
paper, let us study the current in the case of a magnetic 
field that is static and constant inside a finite region and 
vanishes outside of this region. For simplicity we consider 
the magnetic field of an infinitely long solenoid with ra¬ 
dius R that is oriented in the a;3-direction. The magnetic 
field is constant inside the solenoid, points in a;3-direction 
and vanishes outside of the solenoid: 

B [xi,X2,xz) e:iQ{R- (30) 

where 0 is the Heaviside step function, 63 is the unit vec¬ 
tor along the 3rd coordinate axis and B° is the magnetic 
field strength. The Fourier transform of the magnetic 
field configuration (30) can be calculated analytically: 

Bs{k) = B°R^^S{ks) oF, (2-,-^{kl + kl)R^^ , 


( 31 ) 
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FIG. 9: Relative importance of nontrivial corrections to acME (k) due to inter-fermion interactions. Solid lines with symbols 
correspond to the contribution of tree diagram to ocme {k) (first summand in (19)) and dashed lines with the same symbols 
- to the contribution of loop diagrams (second summand in (19)). Since the contribution of loop diagrams turns out to be 
numerically very small, we multiply it by a factor 10^ (rrir = —2.0) or 10^ (m^ = 0.0). 
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FIG. 10: The ratio of the contribution of tree-level diagram 
(first term in (19)) to the contribution of the loop diagrams 
(second term in (19)) to the chiral magnetic conductivity 
o'cME (k) at k — 6 and = 0.05 as a function of the bare 
mass and interaction potential U. 



where qFi is a generalised hypergeometric function [-56]. 

According to the Kubo formula (17) and the expression 
(16) (with an obvious permutation of indices), the Fourier 
transform of the 3rd component of the electric current is 
related to B 3 {k) as 

jsik) = (JCME {k) B^{k). (32) 

Since here we are interested in a rather qualitative es¬ 
timate of the electric current, for simplicity we use the 
continuum expression (28) for acME, which describes our 
numerical data quite well in the phase with unbroken par¬ 
ity (see Fig. 5). 

For experiments it is much more convenient to consider 
the coordinate space current jsix). There is no analytic 
expression for the inverse Fourier transform of (32), and 
we have calculated it numerically. Because of the cylin¬ 
drical symmetry of the magnetic field configuration the 


current jsix) = j^{p) depends only on the radial coor¬ 
dinate p = i^xf + X 2 and on the parameters R, pA and 
B°. The magnetic field strength simply enters the final 
result as a factor, but jslp) is a nontrivial function of pA 
and R. In Fig. 12 we plot jsip) at fixed pA for differ¬ 
ent values of the dimensionless product paR- For small 
PaR the current inside the solenoid is almost constant 

gO 

and its magnitude is close to the value com¬ 

plete agreement with the conventional expression (2). As 
we increase paR the current develops a stronger depen¬ 
dence on p and assumes its maximum near the boundary 
of the solenoid. It is also interesting to note that the 
total current through the {xi,X 2 ) plane is proportional 
to (JCME {k —>■ 0) and hence vanishes. Thus the current 
inside the solenoid is compensated by a “heavy-tailed” 
current density of opposite sign outside of the solenoid. 

To find an estimate for pA in a typical experiment we 
can use the results of [5], where a reasonably good de¬ 
scription of experimental data was obtained under the 
assumption pA ^ p and the chemical potential was 
found to be /X ^ 100 meV. As a reasonable upper 
estimate for the chiral chemical potential we can take 
Pa ~ 10 meV. To observe a strong chiral magnetic cur¬ 
rent inside the sample our numerical calculations sug¬ 
gest that one should have /Xyii? < 1. To achieve this the 
sample size (in the directions orthogonal to the magnetic 
field) should be of order R ~ 10^ eV“^ or, in MKS units, 
R ~ 10“^ mm. This is a reasonable estimate if one takes 
into account the characteristic size of the sample of order 
1 mm and possible inhomogeneities of crystal structure. 

Unfortunately, the characteristic values of the electric 
field strength in the sample are not reported in [-5], and 
hence it is difficult to estimate the stationary value of 
PA directly using the expression (3) derived from the 
anomaly equation. However, it is interesting to note that 
the above discussed dependence of the CME response 
on PaR should result in the nonlinear dependence of 
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FIG. 11: Numerical results for acME (k) at finite temperatures. Point styles denote different sets of model parameters (m^ 
and U) and line styles different values of temperature. The bare chiral chemical potential is ^j!'^ = 0.15 everywhere. 



FIG. 12: Electric current density j'a (p) in the direction of 
the magnetic field as a function of the radial coordinate p for 
different values of the dimensionless quantity paR, where R 
is the radius of the region with nonzero magnetic field. For 
convenience, we show the ratio of js (p) to the conventional 
value j 3 = which follows from (1) and (2). 


the electric current in the sample on the electric field 
strength. Indeed, for sufficiently small electric fields the 
stationary value of pA in (3) and hence the quantity paR 
are small, and the CME response should be well described 
by the conventional expression (2). With the increase of 
the electric field, however, the stationary value of pA and 
hence paR increase, and according to Fig. 12 the CME 
contribution to the conductivity of the sample should de¬ 
crease. 

In our simple estimates we have neglected several com¬ 
plications that can occur in realistic experiments. First 
of all, we have approximated the state of dynamical equi¬ 
librium in the chirality pumping process by a ground 
state of the many-body Hamiltonian, thus completely ne¬ 
glecting all dynamical real-time processes. We have also 
not taken into account that the magnetic field inside a 


sample depends on numerous different parameters and is 
in general a non-trivial function of the applied external 
field. Moreover we did not consider the effect of grain 
boundaries inside the sample or possible contributions 
of boundary states. Nevertheless our simple calculation 
shows that it is possible to observe the chiral magnetic ef¬ 
fect even in the case of static and homogeneous magnetic 
fields and allows us to make educated order-of-magnitude 
estimates that could help in devising new experiments. 


6. CONCLUSIONS 

In this paper we have presented the results of a 
mean field study of a lattice model for a parity-breaking 
Weyl semimetal with on-site inter-electron interactions, 
in which the chiral symmetry is only an emergent low- 
energy symmetry. We have studied the phase diagram of 
the model in the parameter space of bare mass, interac¬ 
tion strength and bare chiral chemical potential. More¬ 
over, for all points of the phase diagram which we con¬ 
sidered we have also calculated the static chiral magnetic 
conductivity within the linear response theory. 

We have found that our model exhibits a phase of spon¬ 
taneously broken CV symmetry (Aoki phase, Axionic in¬ 
sulator phase). The phase transition to this phase is 
of the second order and persists at finite bare chemi¬ 
cal potential and finite temperature. Increasing the bare 
chiral chemical potential leads to a slight decrease in 
the critical interaction strength for the phase transition. 
Our calculations show a strong multiplicative renormal¬ 
ization of the chiral chemical potential in both phases. 
The slope of the renormalized chiral chemical potential 
changes at the phase boundary, but in both phases it 
strongly grows with the interaction strength. This is con¬ 
sistent with a previous study of a model with continuum 
Dirac fermions [27]. 

It turned out that in our model the enhancement of the 
chiral magnetic conductivity (as compared to the non- 
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interacting case) is mostly caused by the growth of the 
renormalized chiral chemical potential with interac¬ 
tion strength. Our estimates of ucme never exceed (by 
absolute value) the universal value with the renor¬ 

malized value of the chiral chemical potential and the 
appropriate number Nj of Dirac cones in the energy spec¬ 
trum. If we adhere to the quantum field-theoretical view¬ 
point on the renormalization, in which only the renormal¬ 
ized quantities are physically observable, we are tempted 
to conclude that inter-electron interactions in fact do not 
result in enhancement of the chiral magnetic conductiv¬ 
ity. To some extent the renormalized value of can be 
indeed observable, e.g., in ARPES studies of the fermion 
dispersion relation. On the other hand, in condensed 
matter systems we can also have an independent way 
of estimating the bare value of fiA, e.g. by using the 
relation similar to (3) in the case of chirality pumping 
process, or from the strain-dependent parameters of the 
tight-binding model [57]. It is therefore a particularly 
interesting question what is the counterpart of the mul¬ 
tiplicative renormalization of fj,A in tbe situation when 
chirality imbalance is created dynamically, and whether 
it leads to some modification of the naive estimates of the 
CME contribution to the longitudinal electric conductiv¬ 
ity based on (3) and (2). 

Furthermore, we have found that spontaneous break¬ 
ing of CV symmetry in the axionic insulator phase (or 
Aoki phase in lattice QCD terminology) results in an 
immediate strong suppression of ucme- This conclu¬ 
sion is in sharp contrast to the behavior of ocme found 
for continuum Dirac fermions with contact interactions 
in [27], where it was found that the CME response is 
particularly large in a phase with broken chiral symme¬ 
try. Thus it seems that the continuum chiral symmetry 
is important for such strong enhancement. It might be 
interesting therefore to consider the CME response also 
for lattice fermions with a continuum chiral symmetry, 
e.g. for staggered fermions [58]. However, in this case 
there are Nf = 4 fermionic species and only one U (I) 
symmetry, which invalidates the direct comparison with 
the analysis of [27]. Yet another possibility would be to 
use overlap Hamiltonian of [59] with Nf = 1 and an ex¬ 
act 17 (I) chiral symmetry. However, from the point of 
view of condensed matter physics the non-local overlap 
Hamiltonian is a rather artificial construction. Moreover, 
chiral symmetry for such a Hamiltonian is generated by 
non-local transformations [59] and is therefore explicitly 
broken by local on-site interactions. 

We have also studied the importance of inter-fermion 
interactions and found that in those regions of the phase 
diagram where the size of the gap in the energy spec¬ 
trum is small, the only practically important effect of 
interactions is the renormalization of the chiral chemi¬ 
cal potential which enters the universal formula (2). The 
relative magnitude of loop corrections does not exceed a 
few tenth of percent along the lines in the phase diagram 
where the gap size is zero. 

We have found that corrections to (Jcme due to inter¬ 


electron interactions are only important when the gap 
size is quite large, and become comparable with the non¬ 
interacting result inside the axionic insulator phase with 
broken CV symmetry (see Fig. 10). Taking into account 
that spontaneous CV symmetry breaking is in fact a re¬ 
mainder of spontaneous chiral symmetry breaking in lat¬ 
tice models with emergent low-energy chiral symmetry, 
this hnding supports the theoretical expectation that the 
relation between acM e and the axial anomaly coefficient 
might break upon chiral symmetry breaking [27, 60]. 
However, it seems that this observation is not very rele¬ 
vant for practical experiments with Weyl/Dirac semimet¬ 
als, since the CME response is very strongly suppressed 
exactly when the corrections to the chiral magnetic con¬ 
ductivity due to interactions become important. 

In addition, we have briefly considered the effects of 
finite temperature on the CME response. In general, we 
found that finite temperature tends to destroy the “axion 
condensate” in the Aoki phase and suppresses the chiral 
magnetic conductivity. 

This work should be considered as the rather crude at¬ 
tempt at the description of the recent experiments with 
Dirac semimetals [4, 5, 61], where chirality imbalance 
is created by chirality pumping in parallel electric and 
magnetic fields, and the CME response results in the 
negative magnetoresistance phenomenon. Since chiral¬ 
ity pumping is a dynamical process, to describe such an 
experimental situation in a more systematical way, one 
should use the real-time linear response theory and adjust 
the mean-held approximation correspondingly. More¬ 
over, one should take into account possible boundary ef¬ 
fects, which are absent when using periodic boundary 
conditions. We leave these technically more advanced 
calculations for further work. 

Finally, let us discuss the assumption of spatially and 
temporally homogeneous Hubbard-Stratonovich helds 
that do not break rotational invariance, which we have 
made in our mean-held analysis. Some theoretical con¬ 
siderations [31, 32, 62] suggest that the system of chirally 
imbalanced fermions is unstable towards the decay of chi¬ 
rality imbalance at the expense of creating helical mag¬ 
netic held conhgurations. In the static limit, which is jus- 
tihable at sufficiently late times after the decay, such spa¬ 
tially inhomogeneous conhgurations should be described 
by the Beltrami equation [33]. In such conhgurations, 
the CME current is expected to vanish [33]. 

Since fermionic condensates in (12) to some ex¬ 
tent can imitate external gauge helds (if, e.g., 
~ OiAi), one can expect that if such in¬ 
stability is indeed relevant, there can exist spatially inho¬ 
mogeneous condensates with lower free energy than that 
of the homogeneous condensates considered in this work. 
Calculations in the holographic Sakai-Sugimoto model 
partially support the possibility of a non-homogeneous 
ground state in the presence of chiral chemical potential 
[35-37], although typically only for sufficiently large val¬ 
ues of fiA of order of the vector meson masses. To check 
for possible unstable directions in the space of all possi- 
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ble configurations of inhomogeneous condensates we have 
investigated the eigenvalues of the Hessian of the free en- 
0$ B 0<i> A ' found no evidence for non-positive 
eigenvalues and conclude that the homogeneous field 
is at least a local minimum in the space of all possible 
fields. Moreover, in a recent work [38] co-authored by 
one of us the real-time decay of chirality imbalance in 
the Hamiltonian similar to (4) was considered, and no 
evidence for the formation of nontrivial spatially inho¬ 
mogeneous gauge field configurations was found. Fur¬ 
thermore, in this work it was found that the decay of 
chirality imbalance is strongly suppressed for small val¬ 
ues of the Fermi velocity. For all these reasons we believe 
that the inhomogeneous ground state is irrelevant for chi- 
rally imbalanced Dirac or Weyl semimetals. 
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Appendix A: Second derivatives of the fermionic 
free energy over the gauge vector potential and the 
Hubbard-Stratonovich field 

In order to calculate the chiral magnetic conductivity 
acME (k) from the Kubo formula (17), we need to cal¬ 
culate second derivatives of the free energy with respect 
to external fields according to the expression (22). Since 
we are working with a finite lattice, we can perform sum¬ 
mation over eigenstates in this expression numerically in 
a straightforward way. 

To this end we need an explicit expression for the 
effective single-particle Hamiltonian. Substituting the 
saddle-point value of the Hubbard-Stratonovich field (23) 
into the definition of the effective Hamiltonian (14), we 
obtain: 


where rjo- is the eigenstate of the operator (Tisin(fci) with 
eigenvalue aS. The corresponding wavefunctions are 
given by: 


4/r(fc) = VsAk) 


Ak-i 




(A4) 


To proceed further, we rewrite (22) in more details: 


o2 jr _ o2 u 

^ = E (•.‘’•Pi ^ I*.‘’.Pi - 

P,S,(T 


-E 




AT cosh' 


(=p#) 




P,S,(7 

q,s' ,(t' 


m If \s,T]^ 


, (A5) 


where the sum is taken over all momenta in the Brillouin 
zone and 




a:,!/ 


Now let us turn to the explicit form of the matrix ele¬ 
ments in (A5). a simple calculation gives: 


f)h 


A=0 


-{.P+5z+iA^’y + 


z+i,y 


), (A7) 


f)h 

^x,y 

d^z,A 


= Sz,J-,yrA, (A8) 


and 


3 

^x,y — ^ ^ i^k^kxy H ^xy 
k^l 


^^x5x,y + iQrUr + (Al) 

and in the absence of external fields in momentum space: 


3 3 

h{k) = E at sin (fci) -I- 270 {hA) + 

i-l i=l 

+7omr -I- Aolsmi + 75Ma- (A2) 

The energy levels of the Hamiltonian (A2) are given by 
formula (25) and its eigenstates have the following form: 


Ts,a{k) = 


les 


W Af-Trir — imi 


<yS-y.A „ 


d^h. 




dAz^-^ j^dAz^^ j 


— ,22 ^ 




P-Sz^.xS 

d^h 
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-PA 
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A=0 
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'f'X.V 




= 0, (AlO) 


= 0,(A11) 


where we denote 


P± 


(Aj ± ^7o 
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(A12) 


Using expression (A6), we can represent the matrix ele- 
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merits in (A5) as: 


dh 




= <Ps,a{p)ji‘Ps',a'i^- 

dh 


3-*(p-9)'2 




d^h 


(s,cr,_^ 




|s,cr,_^ = 


where 


/ ->v / -«i *^21,22 


^ - P+e-^P-^'), 


(A13) 

(AM) 

(A15) 

(A16) 

(A17) 


and e) is a unit lattice vector in the direction i. 

Considering for instance the derivative d’^PojdAidAj, 
we substitute expressions (A13) and (A15) into (A5) and 
perform a Fourier transform with respect to zi and Z 2 ' 

■(«) = 


dAi dAj 


= — y 

r.3 


,ifc(zl-Z2) 


A! 


dAi dA. 


izi,z2), 


(A18) 


and obtain the final expression, which can be easily cal¬ 
culated numerically: 


^3 ^^-^0 


= EE 


dj 


‘Ps,ai]^-^<Ps,aip)n{Ss,a{p)) “ 


rt / 2 \ “r 


^k,0 


AT cosh' 






Ts,aiAjiTs',cT'{q)‘Ps',<7'{^jjTs,a{p) , 

X A {Ss,a{p)) 


(A19) 


ing of the model parameters and observables, which is a 
general feature of the instantaneous potential approxi¬ 
mation for inter-fermion interactions. Therefore in our 
calculations we use vp = A. In this Appendix we provide 
explicit expressions which can be used to restore the de¬ 
pendence on vp in all our results. 

It is easy to check that if we substitute the rescaled 
values 

/vp, /vp, U = U/vp (Bl) 

into the Hamiltonian (4) with the unit Fermi velocity 
^^ = 11^ = 1, we obtain the same Hamiltonian with a 
non-unit value of uf, but rescaled by 1 /uf: 

H(vp = l,A°\pf,u) = 

= —H(vp,m<•°\^if,u)■ (B2) 

Vp \ / 

It is also easy to check that the same rescaling works also 
upon the Hubbard-Stratonovich transformation, if we 
rescale the Hubbard-Stratonovich field as $3, = ^x/vp. 
In particular this implies that if the renormalized model 
parameters /i^, ffir and are known atwF=f’F = lj 
their values at non-unit Fermi velocity can be found sim¬ 
ply as 

mr = vpffir, rrii = vpfhi, ,^A = vpp,A- (B3) 

Substituting the scaling law (B2) into the definition 
P = —Tin ^Tr exp of the free energy, it is 

also easy to obtain 

T(uF,m(°),y,17,T) = 

= vpP(vp = l,A°\llf,U,f), (B4) 

where T = T/vp is the rescaled temperature. Current- 
current correlators which enter the Kubo formula (17) 
are obtained as variations of the free energy over the 
external vector potential A (a;), which is replaced by the 
link phases Ax^k upon Peierls substitution. From the 
explicit form of the Wilson-Dirac Hamiltonian (7) it is 
obvious that Ax^k should not be rescaled when the Fermi 
velocity wf is replaced by uf = 1- We then immediately 
obtain 


where q = p P k and the sum over momentum p is taken 
over the Brillouin zone. 

Final expressions for the derivatives d^P^/dAid^A 
and d^Po/d^Ad^B can be derived in the same manner. 

Appendix B: The role of the Fermi velocity 

In this Appendix we demonstrate that the Fermi veloc¬ 
ity uf in the Hamiltonian (4) amounts to a simple rescal- 


{jx,kjv,i) {vp,m^°\fi‘'2^,U,T'^ = 

5^p(vp = l,m(o),y,i7,T) 

= Vp -^- - — 

= vp{jx,kjv,i) {^F = ■ (B5) 
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Since the coordinates x and hence also the wave vector k 
in the Kubo formula (17) do not scale with Fermi velocity, 
the equation (B5) immediately translates into 

(^CME (k,VF,rn^°\k‘A^= 

= VF crcME {k,VF = ,f^ ■ (B 6 ) 

Using (B 6 ) and (B3), (Bl) one can now eas¬ 
ily see that the linear term in the expansion of 

o'CME {k,VF,'m!''^\ in powers of either the 

bare or the renormalized chiral chemical potential does 
not depend on Fermi velocity. It could be expected, since 
the derivative of acM e over ha is related to the anomaly 
coefficient which also does not depend on the Fermi ve¬ 
locity. 


Appendix C: The fate of the Aoki fingers 


For the sake of completeness we discuss the finite vol¬ 
ume dependence of the “Aoki fingers” for =0. In 
the paper [41] Aoki presents a conjecture for a phase di¬ 
agram of lattice QCD with Wilson fermions. Based on 
calculations in the two-dimensional Gross-Neveu model 
and an effective model of lattice QCD he suggests a phase 
diagram with two phases. In the notation of our paper 
the order parameter for the phase transition is the CV 
breaking mass term mi, which is zero in one phase and 
assumes a finite value in the other phase (Aoki phase, 
Axionic insulator phase). We work with the Hamiltonian 
formalism, where time is continuous and not discretized. 
In the conjectured phase diagram in D dimensions the 
Aoki phase then forms D fingers that touch the bare 
mass axis on characteristic points as the coupling U goes 
to zero. 

We now give a heuristic argument why, in general, the 
fingers do not extend all the way down to touch the bare 
mass axis. At the border of the Aoki phase the pion mass, 
which is given by mj = , has to vanish. In D = d+1 

dimensions the second derivative of the mean-field free 
energy (24) reads 


dmf 


1 1 


E 


y£_i^„(fc) 


_m]_\ 

£_i,„(fc)3y 



where fc is a d-dimensional vector and we sum over all 
k in the Brilloin zone. Let us now consider the limit 
Ls —>■ oo. The sum in (Cl) then becomes an integral: 

(2 ^(“) 

k 

In the case of vanishing Ha'^ the energy levels become de¬ 
generate and we define £_i(fc) := £_i__|_i(fc) = £_i _i(fc). 
Near the phase boundary m^ <C I and around points 


where W -I- rrir = 0 the integral is dominated by the con¬ 
tributions from fc <C I and we can write 


-7 





_ 2 

mf ^ 

{k'^ + j ’ 


where we collect the factors of I/(27r) and the contribu¬ 
tion from the solid angle integral in the (positive) con¬ 
stant 7 . The exact numerical value of 7 is irrelevant for 
the following argument. 
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FIG. 13: Volume dependence of the mean-field phase diagram. 
The Aoki fingers are volume dependent and seem to vanish 
as Ls —>■ 00 . This is only relevant if ~ 0- At finite 
the fingers are not present and phase diagrams for different 
lattice volumes lie on top of each other. In the figure the 
phase boundaries for = 0-10 where shifted by —0.8 in U 
for better visibility. 


As we approach the bare mass axis the term 2/U di¬ 
verges. In order to get = 0 this diver¬ 

gence has to be canceled by the integral. In D = I -|- I 
dimensions the first term diverges like log(mi) and the 
second term becomes simply a constant as mi —^ 0. It 
is possible to cancel the divergence and the Aoki fingers 
can reach all the way down to the bare mass axis. We 
performed numerical calculations in D = 1-1-1 dimen¬ 
sions in the limit —>■ cx) and found that the results are 
consistent with that statement that the fingers touch the 
bare mass axis. 

In D = 3 -I- 1 dimensions, however, the first integral is 
no longer divergent as rrii —> 0. The contribution from 
the second integral is always non-negative and can not 
cancel the divergence. It is therefore not possible to find 
roots of equation (Cl) and the Aoki phase can not extend 
down to the bare mass axis. 

With the same argument it is now straight forward to 
see why the Aoki fingers are volume dependent. If is 
small the sum in (Cl) diverges as IjVk^ for vanishing 
rrii and we can find a solution for jdrnf = 0 even for 
small U. As we increase Lg the sum becomes a better 
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and better approximation of a (non-diverging) integral and (Cl) does no longer have roots for C7 <C 1. 
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